let
  c = [1.,-1.,1.];
  G1 = [-7. -11. -11. 3.;
         7. -18. -18. 8.;
        -2.  -8.  -8. 1.]';
  G2 = [-21. -11.   0. -11.  10.   8.   0.   8. 5.;
          0.  10.  16.  10. -10. -10.  16. -10. 3.;
         -5.   2. -17.   2.  -6.   8. -17.  -7. 6.]';
  h1 = [33. -9.; -9. 26.];
  h2 = [14. 9. 40.; 9. 91. 10.; 40. 10. 15.];
  sol = CVXOPT.sdp(c, [], [], (G1,G2), (h1,h2));
end
